Static overscreening and nonlinear response in the Hubbard Model 



O 

o 

(N 



Erik Koch 

Max-Planck-Institut fur Festkorperforschung, Heisenbergstrafle 1, 70569 Stuttgart, Germany 

(February 1, 2008) 

We investigate the static charge response for the Hubbard model. Using the Slave-Boson method 
in the saddle-point approximation we calculate the charge susceptibility. We find that RPA works 
quite well close to half-filling, breaking, of course, down close to the Mott transition. Away from 
half filling RPA is much less reliable: Already for very small values of the Hubbard interaction 
U, the linear response becomes much more efficient than RPA, eventually leading to overscreening 
already beyond quite moderate values of U. To understand this behavior we give a simple argument, 
which implies that the response to an external perturbation at large U should actually be strongly 
non-linear. This prediction is confirmed by the results of exact diagonalization. 
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I. INTRODUCTION 



Motivated by the surprising accuracy of the random 
phase approximation (RPA) for a half-filled, generalized 
Hubbard model with orbital degeneracy,EJ we ask how well 
RPA reproduces the screening in correlated systems in 
general. While it properly describes the screening when 
the kinetic energy is much larger than the interaction en- 
ergy, RPA becomes wrong in the opposite limit. Its main 
deficiency is that it completely misses the break-down of 
the screening at the Mott transition. Nevertheless, as 
the quantum Monte Carlo calculations in Ref. |] have 
shown, RPA gives a surprisingly accurate description of 
the static charge response on the metallic side of the Mott 
transition, until the system is quite close to the transi- 
tion point. Away from half filling (or, more generally, 
integer filling) one would expect the system to become 
more metallic, and hence, RPA to work even better. We 
will see, however, that this not the case. 

In the present study we use a Slave-Boson methodBi 
to calculate the linear charge response. The advantage 
over more elaborate methods like quantum Monte Carlo 
(QMC) or exact diagonalization are (i) the possibility to 
treat very large systems, such that finite size effects be- 
come negligible, (ii) the efficiency of the method, which 
allows us to thoroughly study the whole parameter space 
of different densities and interactions, and (iii) the pos- 
sibility to directly calculate the linear response, which 
in Lanczos or QMC calculations has to be extrapolated 
from several calculations for different perturbations of fi- 
nite strength. 

There have been already a number of works deter- 
mining the linear response from a one-loop expansion 
around the Slave-Boson saddle-point solutionJaH but 
they mainly focussed on the structure factors and the 
spin response. In contrast, here we are interested in 
the static density response, which can be calculated at 
the saddle-point level, avoiding the problems encoun- 
tered when-sonsidering Gaussian fluctuations of the Slave 
Bosons Jl3c3 In that sense our approach is related to that 



of Refs. [TJI and [jlj, with the main difference that we 
consider the linearized saddle-point equations, in order 
to directly obtain the linear response. As for the accu- 
racy o£-th£-jresults, based on the experience of previous 
worksalljilJ we expect that the static response should be 
well described. 

The model we consider here is the one-band Hubbard 
model with nearest and next-nearest neighbor hopping 
on a square lattice. It is introduced in Sec. II. A. We are 
in particular interested in the response to a point charge. 
In Sec. II. B we describe the approach for calculating the 
charge response using Slave-Bosons at the saddle-point 
level. To check the method and its accuracy, we com- 
pare with the result of exact diagonalization. In Sec. 
II. C we give the results of our calculations: the screening 
as a function of filling and interaction. The most strik- 
ing result is that, already for quite moderate interaction, 
the systems shows a response that is stronger than the 
perturbation (overscreening). The response, of course, 
strongly depends on the doping and changes with the 
next-nearest neighbor hopping, being enhanced close to 
the van Hove singularity. Comparing with the random 
phase approximation, we find that, contrary to our ex- 
pectation, RPA works best close to half filling. In Sec. Ill 
we give an interpretation of these results. Both, doping 
dependence and overscreening can be understood using 
simple large- U arguments. We also discuss why, by the 
same arguments, we do not obtain overscreening in the 
RPA and that overscreening does not mean that the sys- 
tems becomes unstable. The argument for explaining 
the overscreening does, however, imply that the charge 
response of the non-half-filled Hubbard model should be 
strongly nonlinear, which we numerically confirm using 
exact diagonalization. Indeed, we find that the stronger 
the response, the more nonlinear it will be. A conclusion, 
Sec. IV, closes the presentation. 

We finally would like to point out that overscreening 
in the Hubbard model has been observed before and a 
proof, although in quite a different spirit from ours, for 
the necessity, of overscreening for large interactions has 
been given .113 
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II. SLAVE-BOSON CALCULATIONS 
A. Model 

We consider the one-band Hubbard model with near- 
est and next-nearest neighbor hopping (t and t' , respec- 
tively) 



// 'E flafi,a-t' E flk 



U 2J n» " 



on a square lattice with lattice constant a. For the non- 
interacting system (U = 0) the dispersion relation is 



-2t(cos(k x a) + cos(k y a)) — At' cos(k x a) cos(k y a) 



and the density of states has a logarithmic van Hove sin- 
gularity at At'. 

We ask for the response to the perturbation by an ex- 
ternal point charge c at site r c , corresponding to an extra 
term c U n c , CT in the Hamiltonian. Expanding the ex- 
ternal potential in plane waves, we find, given the static 
susceptibility \ q for wave vector q, for the induced charge 
at site Tj 



SUi = 2J Xq SV q 



cU 



E 

q^O 



Xq e 



iq(ri—r c ) 



(i) 



where we have excluded the q = component, since we 
are really interested in the response for a system with 
fixed number of electrons. Below we will particularly 
focus on the response at the site with the point charge 
An = 5rii =c . 

For the noninteracting system the susceptibility (per 

spin) Xq is given by (we give a more general definition, 
which we will need below) 



X, 
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for n — 0. Here N is the number of lattice sites, the 
sum is over all states (first Brillouin zone) and //. = 
(1 + exp(— P(ek — is the Fermi-Dirac function. In 

case of degeneracy the quotient in (||) becomes the dif- 
ferential. In RPA the susceptibility is 



RPA 

A. q 



2 X q 0) 
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B. Method 



As an efficient way for calculating the static response 
we choose the Slave-Boson method, which will allow us 
to easily explore the parameter space (fillings n CT , inter- 
actions U, hoppings t') for quite big lattices (100 x 100) 
and low temperatures (fl — 100/t). We use the formal- 
ism developed by Kotliar and Ruckensteiru (since we are 
only interested in the density response of the paramag- 
net, we need not explicitly take care of the spin-rotation 
invariance) .□ Instead of studying Gaussian fluctuations of 
the Slave- Boson action around the saddle-point (one-loop 
expansion)o~El, we stay at the mean-field level, simply 
linearizing the saddle-point equations for the perturbed 
system. In that sense, our approach is closer to that of 
Rcfs. [l^ and |l4], where the full saddle-point equations for 
a system with finite perturbation were solved. We thus 
avoid the problems- rj&lated to the representation of the 
hopping operator The main advantage of our lin- 
earized saddle-point equations is, that for a plane-wave 
perturbation the resulting linear system becomes block 
diagonal, with identical blocks of size 5x5, independent 
of system size. We expect our method to give a good de- 
scription of the true screening, giv|enJjh&.good experience, 
in particular for static quantitiesH'Ejtll 



We now give a brief outline of the linearization procedure. We start from the Kotliar-Ruckenstein formulation of 
the Slave-Boson methodfl It maps the physical Fermions onto composite particles of pseudofermions and four Bosons, 
representing empty (e,), singly (pj CT ), and doubly (di) occupied sites. Constraints, ensuring consistency between pseud- 



ofermions and Bosons, enforced by Lagrange parameters (A^ and A^), are introduced to eliminate the unphysical 
states in the enlarged Hilbert space. In this enlarged space, the Hamiltonian can be represented in a form, having 
kinetic terms with mixed fermionic-bosonic character (renormalized hopping), while the interaction becomes purely 
bosonic and bilinear. The pseudofermions, merely occurring in bilinear terms, can then be integrated out exactly, and 
one is left with a purely bosonic action integral (including the Lagrange parameters that ensure the coupling to the 
pseudofermions). The simplest approach to the bosonic action integral is the saddle-point approximation, in which 
all Bose fields are time independent. The Slave-Boson action then takes the form 



k 0h 
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where 



1— a 
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and the e nM are the eigenvalues of the renormalized Hamiltonian 



(3) 



The mean-field values of the Bosons and the Lagrange parameters are determined from the saddle-point equations 



r 1 
/r 1 



dS 

dpia 

1 dS_ 

ddi 

dS 

as 



= 2 




<9A 



(2) 





1 


9<W 


1+e 




' ae< 




1 




1+e 




' 9p;„ 




1 


»W 


1+e 




' ddi 




1 




1+e 







For finite temperature /3, this nonlinear system of equations has to be solved self-consistcntly, since the eigenvalues 
e„. CT of the renormalized Hamiltonian (^) depend on the mean-field values of the bosonic fields. The Fermionic terms 
are most easily evaluated in the form 



9/ de /(e) Tr I G a i (e 



c>H CT , 
<9x CT 



where G^e) = (e — Hc)^ 1 is the Greens function. For the homogeneous system, the Boson variables are independent 
of the site index, and solving the saddle point equations for the paramagnet (p = m — p\) essentially reduces to 
self-consistently solving a 3rd-order polynomial in d 2 . 

For calculating the linear response, it is easiest to consider external potentials Vi = SV q cos(qri), which give rise to 



a response: rij 



, + 5n cos(q r*j) and e 2 — e 2 + 2e <5e cos(g r,), 



Expanding to first order and using the relation 



-9 / de/(e)e™E G ^( £ ) cos (l r j) Gj,i( e ) = X 

~ J —OG 



SBMF.n 

1 



cos(grj), 



where Xq BMF ' n is calculated as (||), but using the eigenvalues of the renormalized Hamiltonian (||) for the homogeneous 
solution, we can factor out the coordinate dependence; i.e. we find that the linearized the saddle-point equations are 
block diagonal, with identical 5x5 blocks of the form 
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I -c J 



.denote the partial derivatives of z, and, since we have set p = p a (paramagnetic solution), the 
partial derivative with respect to p means the sum of the partial derivatives with respect to p-\ and p±. Also, we have 



e i ^ee j • • 



introduced A 



A (°> 2 ) 



spin and the Xq° the susceptibilities ([2|) for the noninteracting part of the unrenormalized Hamiltonian taken at the 
renormalized temperature z 2 /3 and Fermi energy (/i — \( 2 ')/z 2 . Solving the linear system, the susceptibility is given 



6e , B = 2( X 



(0,1) 



n a )/z, C = Xq°'^ / z 2 , and D = 2zcq, where eo is the energy density per 



by X, 



SBMF 



dnq/dVq = 4(j)dp/dV q + ddd/dVq), and the response to a point charge follows from eqn. ([!]). 



To check the method, and to verify our code, we com- 4x4 Hubbard model. As can be seen from Fig. |l| the 
pare the results of the Slave-Boson calculation to the lin- Slave-Boson results basically reproduce the response, al- 
ear response calculated by exact diagonalization for a though they tend to slightly underestimate the screening. 
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The good agreement does not come as a surprise.! 
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FIG. 1. Comparison of the linear response calculated by 
the Slave-Boson method with exact diagonalization for a 4 x 4 
Hubbard model (t' = 0) with N a = l, 5 electrons of each spin 
(closed shells) . The dotted lines show the leading order of the 
asymptotic expansion of the Slave-Boson result. RPA is given 
for comparison. 

While the explicit expressions for the susceptibility 

X SBMF are 

somewhat messy, the results of an asymptotic 
expansion for large U/t are fairly compact. We have to 
distinguish two cases: For n a < 1/2, to leading order 



SBMF 
Kq 



-2(1 - r V ) 3 (l - 2n CT )xf 



The important conclusion to draw from the asymtotic 
expansion is that for large U, to leading order, the sus- 
ceptibility is independent of the interaction. Therefore, 
as long as the large-C/ value does not happen to vanish 
(as is the case for n a — 0, 1/2, 1), the linear response to 
a point charge dn/dc will, to leading order, grow linearly 
with U — beyond any limit. 



C. Results 

1. Overscreening 

The results of the Slave-Boson mean-field calculations 
for the response to a point charge as a function of filling 
n a and interaction U for different hoppings t' are shown 
in Figure ^. The most striking feature is that only ex- 
actly at half-filling, the response vanishes (at the Mott 
transition), while for any other filling n a ^ 1/2 the in- 
duced charge density on the site with the test charge 
eventually increases linearly with U . The reason is that, 
in the limit of U going to infinity, the susceptibility of 
the non-half-filled system stays finite, only for n a = 1/2 
does Xq vanish (see Fig. ||). This necessarily means that 
for strong enough interaction and n a ^ 1/2 the induced 
charge will be larger than the perturbation. What is 
quite surprising is, that this overscreening already sets 
in for very moderate values of U . As a function of the 
next-nearest neighbor hopping, we find overscreening for 
(see the first contour line in Fig. |J): 



t'/t 


-0.2 -0.3 -0.4 -0.5 


u/t> 
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( Xq 2) - 4(1 - n CT )e )xf - (Xi - 1 + 2(1 - n.) 2 ) 2 



For fj„ > 1/2 we have 



SBMF 
Kq 



-2nl(2n a 



(x ( q 2) 



4n a e ) X q 0) - ( X \/> - 2n 2 ) 



(i) 



A superficial look at the two expressions might suggest 
that, while they arc almost symmetric under the substi- 
tution 7Z. CT (1 — Tiff), they violate particle-hole symme- 
try. This is, however, not the case. For t' = 0, the results 
are particle-hole symmetric, due to the properties of the 
generalized susceptibilities. To give an impression of how 
the asymptotic expansion works for moderate values of 
U, we included it in Fig. [l]. 



We also remark that, according to the comparison with 
the exact results for a small system (Fig. |l|), the Slave- 
Boson results might even somewhat underestimate the 
true response. 

The second striking feature is the strong doping de- 
pendence of the response. It tends to be strongest for 
quarter filling (rj CT = 1/4 or 3/4), while it vanishes when 
approaching the Mott- or band-insulating regions. This 
doping dependence is somewhat modified by the next- 
neighbor hopping t': The screening becomes enhanced 
in the neighborhood of the van Hove singularity (at At 1 ). 
For t'/t = —1/2, where the van Hove singularity is at the 
lower band-edge, this enhancement is particularly pro- 
nounced, while for t' = 0, where the singularity is at 
half filling, it is masked by the Mott transition. We note 
that around half filling the system could become anti- 
ferromagnetic and that close to the van Hove sipgularity 
and for large U there could be ferromagnetismlla, while 
in our calculations we always restrict ourselves to the 
paramagnetic state. 
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100 x 100 lattice, fit = 100, t'/t = 00.00 

-dn/dc 



100 x 100 lattice, (it = 100, t'/t = -0.10 

-dn/dc 





100 x 100 lattice, fit = 100, t'/t = -0.20 

-dn/dc 




100 x 100 lattice, fit = 100, t'/t = -0.40 

-dn/dc 




100 x 100 lattice, (it = 100, t'/t = -0.30 

-dn/dc 




100 x 100 lattice, (it = 100, t'/t = -0.50 

-dn/dc 




FIG. 2. Linear response to the perturbation by a point charge in the Hubbard model as a function of filling n a and 
interaction U for different values of the next-nearest neighbor hopping matrix element t' . The calculations use the Slave-Boson 
method in the mean field approximation and were performed for square lattices of size 100 x 100 at (3 = 100/t. dn is the electron 
density that is induced on a site with an infinitesimal test-charge dc, i.e. on a site with external potential del] . Contour lines 
are drawn for —dn/dc = 1, 1.5, 2, . . ., i.e. the first contour line marks the onset of overscreening. It is clear that, as the van 
Hove singularity (at At') is shifted to lower energies, the response for less than half filling gets increasingly stronger. For t' = 
the singularity is at half filling and moves towards smaller n CT , until it reaches the lower band-edge for t'/t = —1/2. 
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FIG. 3. Sucsceptibility \ c — X/ 9 ^o X-Q calculated by the 
Slave-Boson method in the limit U — > oo. The calculations 
are for the same systems as in Fig. ^. 



2. Comparison with RPA 



The overscreening obviously has important implica- 
tions for the validity of the RPA. Since in the RPA 
the response can at most equal the perturbation (perfect 
screening), it clearly will not work well, wherever there 
is overscreening. From the direct comparison in Fig. |] 
we find, however, that it already starts to fail for much 
smaller values of U. (The calculations are for a Hubbard 
model with non-zero next-nearest neighbor hopping, to 
avoid perfect nesting.) Surprisingly, RPA seems to work 
best close to half filling, but, of course, only well below 
the Mott transition. Note that the contour lines around 
n a = 1/2 in Fig. |J first mark increasing values, but, when 
approaching the Mott transition, turn negative (see also 
Fig. |). When going away from half filling, where one 
expects the correlated system to become more metallic, 
RPA rapidly fails to give a good description of the screen- 
ing. For better visualization, we give a direct comparison 
of RPA and the result of the Slave-Boson calculation in 
the two regimes (close to and far from half filling) in 
figure ||. For half-filling, RPA only somewhat underesti- 
mates the response for small U, but fails to describe the 
eventual break-down of the screening at the Mott transi- 
tion. Overall, for n a = 1/2 the error in dn/dc is < 0.15 
(relative error < 20%) up to U larger than U c /2. This 
is consistent with what was found in quantum Monte 
Carlo for a half-filled system with orbital degeneracy^! 
although it seems that the orbital degenerate system is 
even slightly better described by RPA. Away from half- 
filling, RPA essentially misses to describe the steep rise 
of the response with U and consequently fails to describe 
the screening even for surprisingly small values of U : Al- 
ready for U > 3 the absolute error exceeds 0.15 (relative 




(dn/dc SBMF -dn/dc RPA ) / dn/dc RPA 



FIG. 4. Difference between the response calculated us- 
ing the Slave-Boson method and RPA for the Hubbard model 
(t'/t = —0.1, to avoid perfect nesting) on a 100 x 100 lattice 
at j3 = 100/t. The upper plot gives the absolute difference 
(contour lines at -0.25, -0.2, -0.15, . . . , 0.25), the lower plot 
the relative difference (contour lines at -0.5, -0.4, . . . , 0.5). 
It shows that RPA works best close to half filling, while it 
quickly gets worse away from half filling. It is to be expected 
that the RPA fails when the overscreening sets in. But it 
turns out that even before that the RPA is already fairly bad. 
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FIG. 5. Comparison of the response calculated by the 
Slave-Boson method with RPA at half and quarter filling. 
The calculations are for the same system as in Fig. ^| (in par- 
ticular, t' 7^ 0, so there is no perfect nesting). 
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III. INTERPRETATION 



A. Doping dependence 



It is clear from Fig. g that the response strongly de- 
pends on the electron density n a . Obviously for the com- 
pletely empty (n a — 0) or the completely filled (n a = 1) 
system there is no response at all, while for half filling 
(n a = 1/2) screening breaks down at the Mott transi- 
tion. In-between the response has a maximum. Since 
this maximum, as a function of electron density, is most 
pronounced for large U, we will consider the screening in 
that limit. For concreteness, we consider a system with 
n a < 1/2. Then, for large U, there will only be empty or 
singly occupied sites. In response to an external pertur- 
bation these electrons and holes will rearrange. Clearly, 
if there are no electrons (n CT close to zero) or no holes 
(n a close to one) there is no room for rearrangement, so 
the response will be very weak. In the other extreme, 
for quarter filling, all electrons and holes can participate 
in the screening. More generally, for the screening of a 
plane wave with nonzero wave vector, the same number of 
electrons and holes have to participate in the response, so 
the response is limited by the density of electrons/holes, 
whichever is smaller. For n a > 1/2 the argument is anal- 
ogous, simply replace empty by singly occupied sites and 
singly occupied by doubly occupied sites. We thus ex- 
pect the response to be strongest around n a — 1/4 and 
3/4, which is roughly what is found in the calculations. 
There is of course also the prominent effect of the van 
Hove singularity, which greatly enhances the response 
and accounts for a shift in the position of the screening 
maximum, as well as for the strong asymmetry around 
half-filling. 



B. Overscreening 

In order to understand the overscreening we also con- 
sider the limit of large U. Again, for a system with less 
than half filling (n CT < 1/2) there will be only singly oc- 
cupied and empty sites. If the test charge is positive, the 
electrons will see a repulsive potential cU. Configura- 
tions with an electron on the site with the test charge 
will therefore be avoided (see Fig. |]): Moving the elec- 
tron from that site to some other empty site results in a 
gain in potential energy c U, at no cost in interaction en- 
ergy. Thus vacating the site with the test charge merely 
costs kinetic energy. Since for large U the potential en- 
ergy will eventually dominate, in the large-?7 limit, the 
site with the repulsive test charge will be empty. This 
implies (i) that for a small perturbation c the response 
will be larger than the perturbation (overscreening) and 
(ii) the response is (almost) independent of the pertur- 
bation (An = n). But this means that the response is 
highly nonlinear: while in linear response one expects 
An/c w dn/dc or, equivalently, An oc c, in the limit 



U — > oo we expect that An is constant, independent of 
c! We will further explore this nonlinearity in the next 
section. 

There are two questions that come to mind. First, also 
RPA describes screening as the interplay between poten- 
tial and kinetic energy. So why is there no overscreening 
in RPA? The answer is that in RPA the electron-electron 
interaction is described by a mean-field: An electron 
interacts with the mean electron density. So in going 
through the argument from above, we have to also ac- 
count for the fact that removing an electron from the site 
with the test charge increases the mean electron density 
on the other sites, consequently increasing the mean-field 
interaction energy. Since, for the infinite system, this in- 
crease scales the same way as the gain in potential energy, 
RPA can give at most perfect screening: For a system 
with M sites, neglecting the kinetic energy, we have 



E/U=\{n-AnY + ^ Y 



1 



M — 1 



+ c(n — An), 



which is minimized for An = (M — 1)/Mc. The sec- 
ond question concerns the stability of the system. The 
perturbing charge induces the system to deviate form its 
homogeneous density in order to counteract the perturba- 
tion. But overscreening means that the resulting electron 
density becomes even more inhomogeneous than the per- 
turbation, while naively one would expect that screening 
tends to restore as much as possible a homogeneous den- 
sity. So does overscreening mean that the system prefers 
an inhomogeneous charge distribution, i.e. that it is in- 
stable against formation of, e.g., a charge density wave? 
The answer is, of course, no. Forming an inhomogeneous 
charge density costs at least some kinetic energy. So only 
when this cost is compensated by a gain in potential en- 
ergy, due to the interaction with the perturbing charge, 
will the charge density become inhomogeneous. 




FIG. 6. Screening of a point charge c in the large-t/ limit 
for less than half filling. Given a positive test charge c the sys- 
tem will gain potential energy c U by emptying the site with 
the test charge. Since there are empty sites in the system, 
no double occupancy needs to be created, so there is no cost 
in interaction energy. Thus vacating the site merely costs ki- 
netic energy. Therefore, if U is very large, an empty test-site 
is favorable: An = n, independent of the perturbation c. 

For an attractive perturbation (c < 0) we can make a 
similar large-C/ argument. Now the system will increase 
the electron density on the site with the test charge. For 
small perturbation the system will tend to put a single 
electron on the site (An = 1 — n). But, as soon as |c| is 
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of the order of unity or larger, it will even pay to create 
a double occupancy, since the gain in potential energy 
(cU) exceeds the cost in interaction energy (U); then 
An = 2 — n. Systems with more than half filling can be 
discussed in just the same way. In that case there is com- 
plete screening for c < by putting an double occupancy 
on the site with the test charge. For 1 > c > the test 
site will be singly occupied, while for c > 1 it will even 
pay to empty it at the expense of creating an additional 
double occupancy. 

At half filling, in the large-[7, limit all sites will be 
singly occupied and for |c| < 1 there will be no response 
at all (An = 0), since moving charge from or to the test 
site would involve the creation of a double occupancy, 
costing an energy U. Clearly, this Mott insulator be- 
havior is missed by the RPA. But it is interesting that 
at half-filling the random phase approximation fails for 
the opposite reason than in the doped case. In RPA the 
screening always involves a cost in interaction energy of 
C/An 2 /2, which, for half filling, underestimates the true 
cost U, while for the doped system it is strongly over- 
estimated. Thus at half filling RPA overestimates the 
response, while for the doped system it is severely under- 
estimated. 

Finally, two points about the arguments given above 
might be worth mentioning. First, the complete- 
screening limits are, of course, nothing but a simple con- 
sequence of the Pauli principle. Second, in our arguments 
we have at no point used the spin of the electrons. Thus, 
as long as there are holes (for n a < 1/2) or non doubly 
occupied sites (for (n a > 1/2) in the system, there will be 
overscreening and nonlinear response, even if the system 
is not paramagnetic. 



C. Nonlinear screening 

As we have argued above, the response to a test charge 
will, in the limit U — ► oo, become practically independent 
of the perturbation: No matter how small the test charge 
c > 0, the system (with n a < 1/2) will respond in the 
strongest possible way (complete screening), by vacating 
the test site (n, = 0). For c < and n a < 1/2 the sys- 
tem will first respond by putting a single electron on the 
test site (n$ = 1), while for c < 1 it will even occupy the 
site doubly (n, = 2). The latter, again, is the strongest 
possible response — complete screening. Thus for large 
U the response becomes more or less independent of the 
perturbation, implying a strong nonlinearity. To check 
this prediction, we have calculated the response to test 
charges ranging from c = — 2 . . . + 1 by exact diagonal- 
ization. 

The results, for a Hubbard model ((' = 0) on a 4 x 4 
lattice with N a = 1 and 5 (closed shell systems), are 
shown in Fig. ^. To emphasize the nonlinearity, we plot 
An/c. If the response was linear, this ratio would be con- 
stant (and equal to the derivative dn/dc). The deviations 



from the horizontal thus show the degree of nonlinearity. 
The horizontal line through unity in the plot marks per- 
fect screening (An = — c), above the response is stronger 
than the perturbation (overscreening). The dotted lines 
give the large-£7 limits discussed above. The curves la- 
beled Hi = and rii — 2 give the complete-screening 
limit. 

Clearly, with increasing interaction U the response gets 
stronger, eventually showing overscreening. But at the 
same time also the nonlinearity of the response (the slope 
of the curve when passing through c = 0) increases. This 
is not surprising, since the curves for finite U are con- 
strained by the response in the large-£7 limit. Given the 
complete-screening limit it is easy to see that overscreen- 
ing is impossible for c > n and c < n — 2. But this implies 
that, the stronger the overscreening, the more nonlinear 
the response has to be. 



5 1-5 



5 1-5 




FIG. 7. Response to a point charge c for a 4 x 4 Hub- 
bard model (t' = 0) with N a =l and 5 electrons (closed shell) 
calculated by exact diagonalization. An is the induced elec- 
tron density on the test-site. The full lines give the results for 
U = 4, 6, 8, 10, 12, 14, 16, with the response getting stronger 
with U. For linear response An/c is independent of c, so 
deviations of the curves from the horizontal indicate the non- 
linearity of the response. The horizontal line through unity 
indicates perfect screening, overscreening above the line. The 
limiting curves (n — m)/c are given by the dotted lines. 
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IV. CONCLUSION 

We have calculated the static density response for the 
one-band Hubbard model using Slave-Bosons at the para- 
magnetic saddle-point. For finite doping we find over- 
screening, i.e. the response exceeds the perturbation. 
This can be understood in the limit of large Hubbard 
interaction U, since strong screening leads to a gain in 
potential energy which increases with U, while it merely 
costs kinetic energy, which is practically independent of 
U. While overscreening thus necessarily occurs for large 
U , we found that it actually already starts at surprisingly 
small U , in particular if the screening is enhanced by the 
van Hove singularity. 

The existence of overscreening might be surprising, if 
one assumes the random phase approximation to give a 
reasonable description of the charge response. In RPA 
there can be no overscreening, because the interaction is 
treated in the mean-field approximation, which implies 
that screening always costs interaction energy, which in- 
creases with U . Therefore the maximum response pos- 
sible in RPA is perfect screening, where the response 
just compensates the perturbation. By construction RPA 
works well when the kinetic energy is much larger than 
the interaction energy. Thus it works well for small U, 
while it necessarily fails at the Mott transition, where it 
fails to describe the break-down of the screening. Sur- 
prisingly, when doping the system away from half filling, 
RPA rapidly get worse, this time severely underestimat- 
ing the response. In that sense the doped system does not 
behave more metallic and RPA-like. Instead its charge 
response more or less resembles what one would expect 
for a Wigner crystal, i.e. in the opposite limit, where 
the interaction dominates over the kinetic energy. Also a 
Wigner crystal should show overscreening, when the crys- 
tal is pinned at the site of a point-charge perturbation, 
at the expense of loosing the kinetic energy associated 
with moving the Wigner crystal as a whole. Clearly, the 
response in that case should also show the same charac- 
teristic nonlincarity as we have found for the Hubbard 
model. _ 

There have been arguments against overscreening ,E3 
based on the intuition that a system with a response 
exceeding the external perturbation would resemble an 
active device, in the sense of electrical network theory, 
and should not be possible. More formal analyses reveal, 
however, that|-o,ve£screening in the charge response is in- 
deed possible,Eir~E!l and does not imply an instability. As 
we have pointed out above, overscreening does not imply 
that the system prefers an inhomogeneous charge den- 
sity. It is only in the presence of the perturbation that 
the cost in kinetic energy for forming the inhomogeneous 
electron density is compensated by the gain in potential 
energy from screening the perturbation. 

An interesting consequence of the overscreening is that 
it changes the sign of the effective interaction between ex- 
ternal charges, i.e. in the case of the Hubbard model, in 



the overscreening regime external charges will experience 
a negative effective U. Moreover, we note that the over- 
screening becomes the stronger the closer the van Hove 
singularity is moved to the lower band edge — in tempt- 
ing analogy with thajnerease of the transition tempera- 
ture in the cuprateso There are, however, two important 
observations to be made about this effective interaction. 
First, of course, it only applies to external charges and 
not to the interaction between the electrons in the sys- 
tem. Second, the concept of an effective interaction only 
makes sense, when that interaction is to a good approx- 
imation independent of the charges under consideration, 
i.e. if the response is linear. 

We have, however, seen that the charge response of the 
Hubbard model can become strongly nonlinear. This can 
again be understood in terms of our large-t/ argument. In 
fact, it follows, that the stronger the response, the more 
nonlinear it has to become. In the extreme case the in- 
duced density becomes independent of the perturbation. 
It is therefore even possible to give strict limits on the 
strength of the perturbation up to which overscreening 
is possible. This nonlinearity has to be kept in mind, 
when using linear response functions, and diagrammatic 
expansions for the Hubbard model. 
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